Transport through cavities with tunnel barriers: a semiclassical analysis 
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We study the influence of a tunnel barrier on the quantum transport through a circular cavity. 
Our analysis in terms of classical trajectories shows that the semiclassical approaches developed for 
ballistic transport can be adapted to deal with the case where tunneling is present. Peaks in the 
Fourier transform of the energy-dependent transmission and reflection spectra exhibit a nonmono- 
tonic behaviour as a function of the barrier height in the quantum mechanical numerical calculations. 
Semiclassical analysis provides a simple qualitative explanation of this behaviour, as well as a quanti- 
tative agreement with the exact calculations. The experimental relevance of the classical trajectories 
in mesoscopic and microwave systems is discussed. 



PACS numbers: 72.23.Ad, 03.65.Sq 



I. INTRODUCTION 

Ballistic transport through quantum billiards has been 
extensively studied in recent years due to its relevance 
for quantum chaos and the possibility of physical appli- 
cations. Realizations of ballistic billiards include struc- 
tured two-dimensional electron gases in semiconductor 
heterostructures and, exploiting the analogy be- 

tween quantum and wave mechanics, microwave cavities 
H . Various experiments have been designed to test the- 
oretical ideas on conductance fluctuations @-^| , weak lo- 
calization [p7| [l0|] and the signatures of classical integra- 
bility. The main theoretical tool for making the connec- 
tion between the quantum and classical properties is the 
semiclassical expansion [0|||. This intuitive and pow- 
erful approach has been tested numerically for the trans- 
port through circular cavities [l^-[T^1. In particular, the 
identification of the most relevant trajectories for trans- 
mission and reflection has been accurately demonstrated 
(analogously to the relationship between the density of 
states and periodic orbits of closed systems ]Tl]]). More- 
over, the semiclassical approach has been extended by 
the inclusion of diffraction effects at the entrance and 
exit of the cavities (ll| . 

In this work we further extend the applicability of 
semiclassical methods in open systems to treat the case 
where tunneling takes place. The modification of the 
trace formula in a closed system by the inclusion of a 
potential step has recently been addressed for a circular 
billiard |l(| in the context of ray splitting. There the 
possibility of electrons entering in the region of higher 
potential has to be taken into account. Our work shows 
that a tunnel barrier within a cavity can very simply be 
incorporated in a semiclassical description and changes 
the relative importance of different classical trajectories 
in a non-trivial manner. Our interest in tunneling in- 
side a cavity stems from a fundamental point of view as 
well as from the fact that experiments with a high poten- 
tial barrier within a quantum dot (pacman) have already 



been performed p7| . The inclusion of a dielectric slab 
within a microwave cavity would also result in a barrier. 

Starting from the well-studied circular billiard [p~3| — P~5|| , 
we introduce a thin barrier placed symmetrically between 
the two leads and extending from the edge to the center 
of the circle as shown by the dashed line in fig. [|. The 
barrier height Vb is variable and allows to interpolate be- 
tween a circle and the billiard studied in ref. fll7f . Of 
special interest will be the regime where the electron en- 
ergy is of the order of the barrier height so that tunneling 
becomes relevant. In the following, we therefore refer to 
this billiard as the circular tunneling billiard. 

We will consider phase coherent and ballistic trans- 
port through the cavity which is attached to two hard- 
wall leads of width W . For a fixed energy E of the inci- 
dent particles there exists only a finite number of trans- 
verse modes N, given by the largest integer smaller than 
(E/Eq) 1 / 2 , which contribute to the transport. Here, 
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is the energy of the lowest transverse mode in the leads 
and M is the mass of the particles. 




FIG. 1. The circular tunneling billiard consists of a circular 
billiard of radius R attached to two opposite leads of width 
W and a thin potential barrier of variable height shown as 
dashed line reaching from the edge to the center of the circle. 
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Within the Landauer formalism |l],|| the two-probe 
conductance g through the structure is just proportional 
to the total transmission coefficient T at the Fermi energy 
E F 
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n.m 



(2) 



The transmission amplitude connecting the incoming 
mode m to the outgoing mode n is given by the pro- 
jection |m 



G nm (x', x,E F ) 



(3) 



dy'dy G(x ', y' , x, y, E F )Xn{y')Xm(v) 



of the retarded Green function G(x' ,y' ,x,y, E F ) of the 
structure onto the transverse modes 
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where we have discarded an unimportant phase factor. 
The longitudinal wave vector is given by 
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The expression (|^) has to be evaluated with x in the in- 
coming and x' in the outgoing lead. For the amplitude 
of reflection one finds the corresponding expression 

ih 2 

Turn — 3nrn ~t~ ^ (&m^n) ^ G nm (x 7 X 7 E F }. (7) 

with x and x' in the incoming lead. 

The transmission and reflection amplitudes can be ob- 
tained numerically by means of the recursive Green func- 
tion method |2(^,^lJ . This method uses a discretized ver- 
sion of the cavity and calculates the Green function by 
starting from the exact Green function in one of the leads 
and successively building up the solution by means of the 
Dyson equation. 

Alternatively, a semiclassical approach to transport 
can be developed from the semiclassical path-integral 
form of the Green function leading to a transmission am- 
plitude Ui 
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given as the sum over classical trajectories s between the 
entrance and exit cross sections with incoming and out- 
going angles 9 and 6' such that sin^ = fhir/kW and 



sin#' = nir/kW (m — ±m, 



±n). The reduced ac- 



tion S is the Legendre transform of the action integral 

S(n, fh, E F ) = S(y' Q , y Q , E F ) + -^(my - hy' ), (9) 

where the starting and end points of the trajectory, 
j/o and y' respectively, are determined by the an- 
gle quantization. For billiards S = kL, where L is 
the length of the trajectory. The amplitude is D = 
(Mv\ cos8'\)- 1 \{dy/d9')e\ and jl is the Maslov index. A 
similar expression holds for the reflection amplitude, with 
the difference that now the trajectories start and end in 
the same lead. 

Direct comparison between the semiclassical ampli- 
tudes (or the conductance) and the exact counterparts 
is rather difficult since the expansion (||) includes an in- 
finite number of terms associated with an exponentially 
large number of contributing classical trajectories. How- 
ever, as in the case of the trace formula, the validity of 
the semiclassical approach can be established by identi- 
fying the Fourier components of t nm with classical tra- 
jectories. This has been done in ref. p5[ and we verify 
it in our system since it gives the starting point of our 
analysis. 

In Sec. H we present the numerical calculation of the 
reflection amplitude and its interpretation in terms of 
classical trajectories without and with a high barrier. In 
Sec. Ell we consider the influence of a tunnel barrier on 



the different trajectories contributing to the reflection 
and transmission amplitudes and find a nonmonotonic 
behaviour. This behaviour is modeled in Sec. |y| within 
a modification of the semiclassical transmission ampli- 
tudes that includes tunneling in a very simple way. In 
Sec. we consider a similar analysis for the transmission 
and reflection coefficients, and show that the analysis be- 
comes considerably more involved since we now have to 
deal with pairs of trajectories. In the final section we 
present our conclusions and discuss the extension of our 
work to the case with magnetic field. 



II. IDENTIFICATION OF CLASSICAL 
TRAJECTORIES 

In fig. H we present the (exact) total reflection coef- 
ficient together with the contribution from the lowest 
mode |r*n | 2 . Obviously, for 1 < kW/n < 2 the two 
quantities coincide since there N = 1. In order to iden- 
tify the classical trajectories we carried out a discrete 
Fourier transformation of 600 values of the complex re- 
flection amplitude r\\ calculated over a momentum in- 
terval ranging from 1 to 7 it /W . The ratio between the 
radius R of the circle and the width W of the leads was 
in all calculations taken to be R/W = 3. 
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FIG. 2. Total reflection of the circular billiard and the 
squared modulus of rn as a function of the Fermi momen- 
tum. For 1 < kW/n < 2, both reflection coefficients coincide 
as there is only one propagating mode. 



Pi 




FIG. 3. Length spectrum TZn for (a) the case without a 
barrier and (b) with an infinitely high barrier. Lengths are 
scaled to the radius R of the circle. In the absence of a bar- 
rier we identify the peaks in the length spectrum with periodic 
trajectories. 

The semiclassical form (||) of the transmission ampli- 
tude implies that its Fourier transform with respect to 
the momentum k — (2ME-p/h 2 ) 1 ^ 2 should exhibit peaks 
at lengths L corresponding to the contributing classical 
trajectories. Similarly, the reflection amplitudes can be 
interpreted in terms of classical paths. For a detailed 



analysis we now consider the power spectrum of rn with 
respect to length. This quantity, which is shown in hg. ||a, 
is given by the squared modulus of the Fourier transform 
of rn and will be denoted as TZu(L) in the following. 

In complete agreement with ref. |15| , we can establish 
a correspondence between peaks of IZn and the classical 
trajectories including their repetitions. The first peak is 
not a classical trajectory contributing to reflection, but 
corresponds to diffraction off the lead mouths p3| , p^[ . 
This effect can be interpreted in terms of a trajectory 
that gets reflected back at the right lead (ghost path). 
For larger lengths L, we can identify a triangular path, a 
five-star path, a seven-star path, and so on. In agreement 
with the semiclassical quantization of the initial and fi- 
nal angle, the star-shaped paths are the most important 
ones for small mode numbers which favor the forward 
direction. 

The resolution of the length spectrum is limited by 
the width of the momentum interval used for the Fourier 
transformation. While this can easily be controlled, there 
are also intrinsic effects restricting the resolution. Within 
a semiclassical picture, at given energy and mode num- 
bers the angle quantization selects the appropriate paths. 
Depending on the width of the leads, the transverse po- 
sition of the starting and end points in the leads are vari- 
able and therefore a given type of trajectories exists in 
a certain momentum interval. In the Fourier transfor- 
mation these trajectories will contribute with different 
weights since the corresponding action will depend on 
the momentum, thus yielding a finite resolution. Already 
the fact that a type of trajectory effectively contributes 
only in a finite momentum interval may limit the reso- 
lution more strongly than the finite interval imposed by 
the numerics. In this respect it is also important to note 
that quantum mechanical diffraction effects at the lead 
mouths [jl3 15| influence the effective momentum interval 
and may be relevant for the resolution. 

Placing a sufficiently high barrier into the cavity yields 
a reflection coefficient (not shown here) uncorrelated to 
that presented in fig. 0. On the other hand, we ex- 
pect that individual trajectories would be greatly altered 
by the barrier and thus we analyze the length spectrum 
lZii(L,Vb) as a function of the barrier height. For very 
high barriers one can see in fig. |b that new length scales 
have appeared rendering the identification more involved 
as compared to the case of vanishing barrier. 

Some features are easily explained like the appearance 
of a peak at length 2R, while in the absence of a barrier 
the minimum length is 4i?. In the presence of a high 
barrier, the direct path may get reflected at the barrier 
thus leading to a peak at half of the previous minimal 
length. The peak at AR now consists of two contribu- 
tions, namely the direct path which is reflected at the 
right lead and twice the direct path reflected at the bar- 
rier which involves a reflection at the left lead mouth. 
The hierarchy continues with a smaller peak at about 
6R which stems from three repetitions of the direct path 
reflected at the barrier. 
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In principle, it is not clear that an analysis in terms of 
classical paths is applicable for arbitrary barrier height 
since tunneling necessarily implies non-classical trajecto- 
ries. However, we will show that such an analysis is still 
possible and helpful towards the understanding of the 
transport problem. For instance, comparison between 
figs. [|a and b shows a large suppression of the harmonic 
coming from the triangular path, while the five-star path 
component is much less affected. Simple arguments given 
in the next section explain this difference in behaviour. 



III. PATHS IN THE PRESENCE OF A TUNNEL 
BARRIER 

For a semiclassical analysis of the energy-dependent 
transmission and reflection spectra, we first have to dis- 
cuss how the classical paths are modified by the barrier. 
The cases where well-defined classical paths exist, are 
those of the circular billiard (V& = 0) and the circular 
billiard with a very high barrier (VJ, — oo). While post- 
poning a more detailed analysis to Sec. we expect that 
at intermediate barrier heights, the transmission and re- 
flection amplitudes for the billiard should be given by 
both classes of trajectories properly weighted according 
to the transmission and reflection coefficients of the bar- 
rier. The length spectrum of the reflection amplitude, 
referred to as reflection spectrum in the following, of the 
circular billiard shown in fig. [|a displays distinct peaks 
which can be associated with a triangle, a five-, and a 
seven-star. In the following discussion we will focus on 
these three trajectories. 

We start with the triangular path as the simplest case. 
In the absence of a barrier, the trajectory just follows the 
triangle as shown in fig. [|a. As the barrier height V& is 
increased the transmission probability through the bar- 
rier decreases and for high barriers the original triangle 
is no longer a possible path. Accordingly, the peak in the 
reflection spectrum corresponding to the triangle will de- 
crease in amplitude with increasing barrier height. For 
sufficiently large Vb, the possibility of reflection at the 
barrier has to be taken into account. As a consequence 
of this reflection the path will no longer continue on the 
original triangle shown as dotted line in fig. ^ but follow, 
at least for sufficiently thin barriers, the dashed-dotted 
line obtained as mirror image with respect to a vertical 
line through the barrier. It is important that this path 
has the same length as the original path. However, now 
the end point no longer lies in the entrance lead but in the 
opposite lead and thus the path reflected at the barrier 
will contribute to the transmission through the billiard. 
Correspondingly, the triangle will become more impor- 
tant in the transmission spectrum as the barrier height 
is increased. This qualitative discussion is confirmed by 
the numerical results for the reflection and transmission 
spectra, IZn and 7~u, shown in fig. |B| as diamonds and 
triangles, respectively. 





FIG. 4. The triangular path in the circular tunneling bil- 
liard may follow the dotted triangle or its mirror image shown 
as dashed-dotted line, a) Trajectory (solid line) which is 
transmitted at the barrier and contributes to the reflected 
trajectories in the billiard, b) Trajectory (solid line) which is 
reflected at the barrier and therefore contributes to the trans- 
mission through the billiard. 
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FIG. 5. Variation of the peak height in the length spectrum 
for the triangular path as a function of the barrier height 14. 
The diamonds and triangles correspond to TZn and 7n, re- 
spectively. Both quantities are normalized with respect to 
7?.n(0), i. e. the case of vanishing barrier. 

While the triangular path cannot appear in the reflec- 
tion spectrum for very high barriers, it is interesting to 
note that there is a peak at a length corresponding to 
two repetitions of the triangular path. This can readily 
be verified by comparing figs. |^a and b. While one re- 
flection at the barrier changes the exit lead from the left 
to the right, an additional reflection restores the left lead 
as exit lead. Therefore, the peak can be identified with 
two repetitions of the path shown in fig. ^b including a 
reflection at the right lead due to diffraction at the lead 
mouth. 

The behaviour of the five-star trajectory is more com- 
plex due to the fact that it crosses the barrier twice. 
Making use of the same geometrical arguments as for 
the triangle, we may distinguish four different classes of 
trajectories shown in fig. n which correspond to two trans- 
missions at the barrier (a), one reflection and one trans- 
mission (b and c), and two reflections (d). In fact, since 
the trajectory may either start into the upper or lower 
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FIG. 6. Four different classes of paths related to the 
five-star: paths with two transmissions or two reflections at 
the barrier (a and d) contribute to the reflected trajectories in 
the billiard while paths with one transmission and one reflec- 
tion at the barrier (b and c) contribute to the transmission 
through the billiard. 

half of the billiard, these four classes correspond to eight 
different trajectories which can be obtained by reading 
the diagrams in the forward and backward direction. 

Like for the case of the triangle, a reflection at the 
barrier changes the sense of rotation in the circle thereby 
changing the exit lead. As can be seen from fig. @, the 
paths with an even number of reflections (a and d) con- 
tribute to the reflected paths through the billiard while 
the paths with an odd number of reflections (b and c) 
contribute to the transmission through the billiard. 

Let us first consider the trajectories contributing to 
the reflection. The trajectory (a) will only contribute 
for small barrier heights because it has to be transmit- 
ted through the barrier twice. On the other hand, the 
trajectory (d) will appear only for rather high barriers 
since it requires two reflections at the barrier. As a con- 
sequence, we expect that the peak in the reflection spec- 
trum corresponding to the five-star trajectory will exhibit 
a minimum at intermediate barrier heights where none of 
the trajectories (a) and (d) contribute significantly. The 
other trajectories (b and c) appear in the transmission 
spectrum at intermediate barrier heights because they 
have to be reflected as well as transmitted once at the 
barrier. 

We conclude from the discussion of the five-star trajec- 
tory that in general the dependence on barrier height of 
the peak heights in the transmission and reflection spec- 
tra should be nonmonotonic. This is confirmed by the 
numerical data shown as diamonds {7tn) and triangles 
(Tn) in fig. 0. The lines shown there are results of a 
semiclassical analysis which will be discussed in detail in 
the next section. 

For this nonmonotonic behaviour it is crucial that 



FIG. 7. Variation of the peak height in the length spec- 
trum for the five-star trajectory as a function of the barrier 
height Vb- The diamonds and triangles correspond to TZn and 
lii, respectively. Both quantities are normalized with respect 
to 7?.n(0), i. e. the case of vanishing barrier. The solid and 
the dotted line give the results of a semiclassical calculation 
taking into account tunneling through the barrier. 

the reflection of a classical path at the barrier does not 
change the length of the trajectory. Only then a change 
in the barrier height will not affect the position of the 
peak and contributions of different paths have to be co- 
herently superposed. 

The behaviour of the peak heights becomes more com- 
plex as the trajectories encounter the barrier more often. 
This will become clear from our final example, the seven- 
star trajectory. In this case, the trajectory encounters the 
barrier three times giving rise to eight classes of trajecto- 
ries. We may classify these trajectories according to their 
behaviour at the barrier by assigning a T or an R for each 
transmission or reflection, respectively. Then, the paths 
contributing to the reflection spectrum are those contain- 
ing an even number of R, namely TTT, RRT, RTR, and 
TRR. While the first trajectory will contribute for very 
small barriers, the other trajectories appear only for suf- 
ficiently high barriers. At very high barriers none of these 
paths is allowed. Accordingly, the peak in the reflection 
spectrum associated with the seven-star will initially de- 
crease with increasing barrier height, exhibit a minimum 
followed by a maximum and then go to zero as the barrier 
height becomes very large. This behaviour can readily be 
verified by comparison with the numerical results for TZ\\ 
shown as diamonds in fig. |[ 

The behaviour of the transmission spectrum 7~n shown 
in this figure as triangles can be understood along the 
same line of reasoning. The trajectories contributing to 
this spectrum are those with an odd number of reflections 
at the barrier, i.e. RTT, TRT, TTR, and RRR. Since 
there is at least one reflection at the barrier, seven-star 
trajectories may contribute to the transmission spectrum 
only for finite barrier heights. With increasing height 
there will be a maximum followed by a minimum and at 
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FIG. 8. Same as in fig. R, but for the seven-star trajectory. 




FIG. 9. Trajectory corresponding to a half of an eight-star 
trajectory, a) The trajectory starting into the lower half of 
the billiard contributes to the transmission through the bil- 
liard for arbitrary barrier height, b) The trajectory starting 
into the upper half of the billiard contributes to the transmis- 
sion in the absence of a barrier (left) and to the reflection for 
high barrier (right). 

very high barriers the RRR-trajectory will contribute. 

So far, we have concentrated on trajectories which 
in the absence of a barrier contribute to the reflection. 
The behaviour of classical paths connecting two different 
leads, i.e. paths contributing to the transmission, as a 
function of the barrier height is slightly more complex. 
In this case, one has to distinguish the paths first go- 
ing into the upper and lower half of the billiard. As an 
example, we consider the paths shown in fig. ^ which 
represent one half of an eight-star. From this figure it 
becomes clear that depending on the incident angle of 



the path there will be an even or odd number of contacts 
with the barrier. Accordingly, the path shown in fig. |^a 
will behave very much like the five-star which also en- 
counters the barrier twice. On the other hand, the path 
shown in fig. ^|b has just one contact with the barrier and 
its dependence on the barrier height resembles that of the 
triangle. As a consequence, while both trajectories will 
contribute to the transmission for low barriers, only one 
of them (fig. ^a) will do so at high barriers. The other 
path will turn into a reflected path (right part of fig. ^|b) 
instead. 

This line of reasoning generally applies to trajecto- 
ries which are transmitted in the absence of a barrier. 
Even though the analysis now becomes more compli- 
cated, these paths also lead to nonmonotonic reflection 
and transmission spectra as a function of barrier height. 
Again an increasing number of contacts with the barrier 
will lead to an increasingnumber of extrema. An exam- 
ple will be shown in Sec. [V] where the spectra of the total 
transmission are discussed (fig. |l3| ). 

The previous qualitative discussion allowed us to un- 
derstand the effect of a tunnel barrier by simple consid- 
eration of classical trajectories. In the next section we 
will show that a quantitative agreement with the exact 
calculations can be obtained within a semiclassical ap- 
proach where the possibility of transmission or reflection 
at the barrier is incorporated. 

IV. SEMICLASSICAL DESCRIPTION OF THE 
CIRCULAR TUNNELING BILLIARD 

We now want to include tunneling into the semiclassi- 
cal picture while remaining rather close to the expression 
for the semiclassical Green function (||) in terms of clas- 
sical paths. To this end, we multiply the contribution of 
the classical paths by amplitudes at or a r accounting for 
each transmission or reflection of the classical path at the 
barrier. 

The treatment of a barrier of finite length in the cir- 
cular tunneling billiard represents a rather complicated 
two-dimensional problem. However, we may approx- 
imately describe the behaviour of an electron at the 
barrier as plane wave encountering a barrier of infinite 
length. Then the problem may be separated into the 
directions perpendicular and parallel to the barrier and 
the only parameter describing the scattering geometry is 
the incident angle 4>. At this point it is important to 
note that the sequence of transmissions and reflections 
at the barrier matters. For example the five-star trajec- 
tories shown in figs. |^b and c which, if read from left to 
right, correspond to TR and RT, respectively, have dif- 
ferent incident angles for the transmission and reflection 
events. 

For an infinitely long barrier the relevant momentum 
component is the one perpendicular to the barrier 

k± = kcos(cj)). (10) 



G 



We now may use the standard results for one-dimensional 
barrier penetration to approximate the transmission and 
reflection amplitudes by 



Ot(<f>) = 



2k±k'j_ exp(— ik±b) 



2k_Lk' x cos(k' x b) - i(k 2 j_ + k' 2 ) sin(fc^6) 



k' 2 - k 2 
<Xt{4>) = i 2k_ L k' 1 ~ sin ( fc ^ b ) at ' 



respectively, where 



k' 



2MV h 
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(11) 



(12) 



(13) 



and b is the width of the barrier. Expressions ([ll]) and 
( |l2| ) reduce to the usual one-dimensional expressions for 
incident angle = 0. 

We emphasize that it would not be sufficient to take 
the modulus of at and a r since in general the contribu- 
tions of different paths have to be added up coherently. In 
fact, destructive interference of paths is quite important 
for the interpretation of the barrier height dependence of 
the length spectra. As can be seen from (O) , the relative 
phase factor between the reflected and the transmitted 
path is always ±i. Since changing a reflection at the bar- 
rier into a transmission and vice versa will change the 
exit lead, two classical paths going to the same exit lead 
differ by their behaviour at an even number of barrier en- 
counters. Therefore, the relative phase factor will be ±1. 
In addition, a reflection at the barrier does not change 
the classical amplitude D, so that the two paths either 
interfere perfectly constructive or destructive. 

As an example we consider the minimum in the peak 
height of the reflection spectrum corresponding to the 
five-star shown in fig. [7] There is such a pronounced 
minimum only because the two contributing paths, TT 
and RR (cf. figs. ^|a and d), interfere destructively. On 
the other hand, TZn is not vanishing at the minimum. 
This is due to the fact that the Fourier transformation 
has to be taken over a finite energy interval. Since the 
barrier height at which the minimum occurs is energy- 
dependent, the minimum of TZn will be smeared out. 

For finite barrier width, the phase factors appearing in 
the transmission and reflection amplitudes have an ad- 
ditional effect. The contribution of the barrier region 
to the total action of the path will depend on the bar- 
rier height which will result in an effective change of the 
length of the trajectory. This becomes more important as 
the width and height of the barrier are increased. How- 
ever, for the parameters used here, the barrier is thin 
enough so that the change in length is below the resolu- 
tion of the discrete Fourier transformation. Nevertheless, 
the peak height is affected. This may become important 
for sufficiently high barriers and cause the decrease in 
the reflection spectrum of the five-star (fig. ^) at large 
Vb- For high barriers of finite width slight changes in the 
scattering geometry with respect to the ideal geometry 
for a thin barrier may also affect the peak height. 



We now turn to a more detailed discussion of the quan- 
tum mechanical results and those obtained from the semi- 
classical approach just introduced. The data are com- 
pared in fig. ^ for the five-star trajectory and in fig. |^ for 
the seven-star. In both cases the diamonds and triangles 
correspond to the quantum mechanical results for IZn 
and Tn, respectively, while the solid and dotted lines are 
the corresponding semiclassical results with the modifi- 
cations described above. 

The expressions for the reflection and transmission am- 
plitudes ( |lT| ) and ( p"2| ) depend on the momentum compo- 
nent k± perpendicular to the barrier and the height Vb 
and width b of the barrier. The geometry of the classical 
path together with the Fermi energy determines k±. In 
the quantum calculations the barrier was implemented 
by increasing on three lattice points the potential to Vb. 
For the curves shown in figs, [t] and ^ we used an effective 
barrier width of 3.5 lattice spacings. 

From figs. and ^ we find, that at not too high barriers 
the agreement between the quantum mechanical data and 
the semiclassical theory modified for tunneling is very 
good. On the other hand, for high barriers deviations do 
occur. This seems to be at odds with the fact that in 
the limit of infinite barrier the modified theory becomes 
equivalent to the usual semiclassical expansion. 

A qualitative deviation appears for rather high bar- 
riers in the transmission spectrum of the five-star and 
the reflection spectrum of the seven-star where the quan- 
tum mechanical data saturate at a finite value. This has 
been checked for barriers as high as 703-Eo- On the other 
hand, the semiclassical result decreases to zero since at 
least one barrier transmission is needed in order to get 
a path which describes transmission through the billiard 
and has the length of a five-star. The same holds for 
paths which describe reflection at the billiard and have 
the same length as a seven-star. 

This discrepancy may be explained by paths with 
length close to those of the five- or seven-star. As an 
example we consider the five-star for which the reflection 
spectrum in the absence of a barrier and the transmis- 
sion spectrum for very high barrier are shown in fig. [l^ 
as dashed and solid lines, respectively. The highest peak 
of the dashed line indicates the length of the five-star. As 
expected, there is no peak at this position in the trans- 
mission spectrum. However, there exists a peak at some- 
what smaller length which is broad enough to yield a con- 
tribution at the length of the five-star. The correspond- 
ing classical trajectory is shown in fig. [ll]. It involves 
three reflections at the barrier as well as a reflection at 
the exit lead, a quantum mechanical effect (l5| already 
mentioned above. 

These arguments also provide a partial explanation of 
the quantitative deviations found for high barriers in the 
reflection spectrum of the five-star trajectory and the 
transmission spectrum of the seven-star. In addition, as 
discussed in the previous section the discreteness of the 
Fourier transformation in combination with the change 
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FIG. 10. Length spectra in the vicinity of the five-star 
peak. The dashed line corresponds to Tin in the absence 
of a barrier while the full line corresponds to 71 1 for a high 
barrier. The scales for 7u and TZn differ by a factor of three. 




FIG. 11. Trajectory appearing in the transmission spec- 
trum for high barriers at a length close to that of the five-star 
trajectory. 

in path length as a function of barrier height results in 
another source of discrepancy. 

While the agreement between the quantum mechani- 
cal and the semiclassical calculation demonstrates that 
for prominent peaks like the ones corresponding to the 
five- and seven-star the interpretation in terms of classical 
paths is possible, the deviations just discussed show that 
nevertheless the identification is not necessarily straight- 
forward and requires a certain amount of caution. The 
dependence of peaks in the length spectrum on the bar- 
rier height may be of help in the identification since the 
number of extrema is related to the number of barrier 
encounters of the corresponding classical trajectory. 

V. TOTAL TRANSMISSION AND REFLECTION 
IN THE CIRCULAR TUNNELING BILLIARD 

The scattering amplitudes discussed above should be 
accessible to microwave experiments where the "Fermi 
energy" can be varied by changing the input frequency. 
A barrier could be introduced by placing a dielectric slab 



into the cavity. On the other hand, in transport experi- 
ments on ballistic microstructures it is the conductance 
which is measured. According to the Landauer formula, 
eq. (|^) , the conductance through a cavity is proportional 
to the total transmission coefficient. Recently, measure- 
ments of transport through billiards with variable Fermi 
energy have been performed ]7],^2| so that it appears fea- 
sible to experimentally determine the length spectrum of 
the total transmission for the circular tunneling billiard. 
Since the semiclassical approaches developed so far pro- 
vide a basis for the understanding of quantum transport, 
we would like to analyze now the effect of a barrier on 
the scattering probabilities. 

The analysis for the transmission and reflection prob- 
abilities is more difficult than for the amplitudes because 
taking the squared modulus makes |t nm | 2 depend on pairs 
of trajectories. Accordingly, a Fourier transformation 
yields peaks at lengths which correspond to differences 
between the lengths of two classical paths. Additionally, 
in order to calculate the total reflection or transmission, 
one has to take into account all scattering amplitudes r nm 
or t nm where the mode number is restricted by the Fermi 
energy via the maximum mode number N introduced in 
Sec. § This will increase the number of different paths 
to be considered since at higher mode numbers paths en- 
closing a larger angle with the lead axis become relevant. 

Moreover, when a new mode opens up it starts by being 
completely reflected, and the reflection coefficient jumps 
by one. This staircase effect will introduce high harmon- 
ics in the reflection coefficient which are not related to 
classical trajectories. Therefore we will focus our atten- 
tion on the total transmission coefficient T. 

As an example, the length spectrum of T for the cir- 
cular billiard without barrier is shown in fig. [l^. Due to 
the large amount of possible length differences the spec- 
trum is quite complex. Most of the peaks correspond 
not only to one pair of trajectories but to a combination 
of several pairs of approximately the same length differ- 
ence. For example the highest peak, which is found at 
L = i.lR, contains contributions of four different tra- 
jectories, namely a half of an eight-star (cf. fig. of a 
twelve-star, a sixteen-star, and a twenty-star. From these 
trajectories, three length differences, 4.20i?, 4.10i?, and 
4.06i?, can be constructed which are all quite close to 
L = 4.1-R. In principle, higher-order stars lead to fur- 
ther length differences in the same range. However, their 
contributions are negligible. Since the resolution of our 
discrete Fourier transformation is 0.1R and the typical 
peak width is twice as large (cf. fig. ||), the contributions 
from the different combinations of trajectories cannot be 
resolved. 

Given the complexity of the length spectrum, we may 
ask the question of whether the inclusion of a barrier 
can be treated as a perturbation that simply randomizes 
the phases of the contributing trajectories as do shape 
distortions (|,|2^,|24|] or magnetic field changes (JJlJ] . Our 
analysis in the previous sections suggests that this is not 
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FIG. 12. Length spectrum of the total transmission T for 
a circular billiard without barrier. 
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FIG. 13. Variation in the peak height of the transmission 
spectra T (diamonds), Tn (triangles), and T22 (stars) at 
L = 4.1i? as a function of the barrier height V&. The data are 
normalized to the respective values at Vj, = 0. The joining 
lines are only guides to the eye. 

the case. The inclusion of a barrier has the effect of sup- 
pressing the contribution from certain trajectories while 
increasing that of the symmetry related ones, resulting in 
the nonmonotonic behaviour of the peaks in the length 
spectrum as a function of the barrier height which cannot 
be interpreted as conductance fluctuations. 

In fig. 13 we present the dependence of the length spec- 
trum on the barrier height for the total transmission co- 
efficient as well as two individual probabilities, |tn| 2 and 
|*22 1 2 - In the following, these spectra are referred to as 
T, Tn, and T 2 2, respectively. The data are given for 
the peak at L = 4.1R to which several paths contribute 
as discussed above. 

The nonmonotonic behaviour resulting from the de- 
pendence of the individual trajectories (of the contribut- 
ing pairs) on Vb illustrates the interplay between semi- 
classics and tunneling. The fact that T as well as Tn and 
T22 roughly exhibit the same structure implies that the 
same pairs of paths contribute to various transmission 



probabilities provided the momentum interval is large 
enough. Of course, predicting the peak height depen- 
dence as we did for the scattering amplitudes is more 
difficult since the relative weight and the phases of the 
trajectories of the contributing pairs become relevant. 
We therefore do not attempt to attain the agreement of 

Sec. Irvl. 



Instead, we will give a qualitative explanation for the 
difference between Tn and T22 ■ With growing number of 
reflections at the circle, the incident angle of the half-star 
trajectories decreases, thus suppressing the contribution 
of the higher-order stars, especially to T22- Since the 
higher-order stars lead to more extrema in the Vb depen- 
dence, the number of extrema of T22 should be smaller 
than that of Tn. This is consistent with the numerical 
results of fig. [b3[ The number of extrema suggests that 
the dominant contribution to T22 stems from the combi- 
nation of one half of an eight- and of a twelve-star. One 
half of a sixteen-star should yield a contribution with a 
second minima for which there is no clear indication. On 
the other hand, Tn displays a second broad maximum 
which is probably a combination of all the extrema from 
higher stars. 

Obviously, the experimental resolution is limited and 
it would not be possible to observe a length spectrum 
as structured as that of fig. O. In mesoscopic systems 
temperature plays two roles. On one hand, it controls 
the inelastic processes like electron-electron and electron- 
phonon interactions. This results in a cut-off length given 
by the inelastic mean-free path L$ beyond which no tra- 
jectory should contribute in a semiclassical expansion. 
Therefore, length differences between trajectories which 
are longer than L<$, should not be considered, render- 
ing in practice a length spectrum smoother than that of 
fig. [l^, thus allowing for a simpler analysis. In the lan- 
guage of ref . [jl5| , the inclusion of L$ corresponds to the 
"smoothing-then squaring" process, since it is the scat- 
tering amplitudes which loose their high harmonics. The 
other effect of temperature is the rounding of the Fermi 
surface that, within our semiclassical approach, cuts the 
large length differences. That is, we should not consider 
the high harmonics of fig. 12 beyond the thermal length 
Lt- Microwave cavities are not sensitive to Lt, but have 
a finite L$ , which can reach very large values in state-of- 
the-art experiments [pi. 



VI. CONCLUSIONS 

In this work we have studied quantum mechanically 
and semiclassically the effect of a tunneling barrier on 
quantum transport through ballistic cavities. This ef- 
fect is most evident in the length spectra, i.e. the Fourier 
transform of the energy-dependent transmission and re- 
flection amplitudes of the cavity. We have shown that 
the peak heights in the quantum mechanical length spec- 
trum exhibit a nonmonotonic variation upon increasing 







the barrier height. This behaviour is quantitatively re- 
produced by combining a semiclassical approach to con- 
ductance with a simple model for barrier reflection and 
transmission of paths. The model provides furthermore 
an intuitive physical picture of the underlying process 
leading to the variations in the peak heights: They reflect 
the superposition of coherent contributions from paths 
being reflected at or transmitted through the barrier. 

This mechanism is a clear-cut manifestation of tunnel- 
ing orbits in the conductance of quantum billiards, which 
should be in principle observable in experiment. The ef- 
fect of the tunnel barrier on individual peaks in the length 
spectrum manifests itself in the energy-dependent trans- 
mission and reflection coefficients. 

We have also performed a corresponding analysis of 
the effect of a tunnel barrier on the area spectrum, the 
Fourier transform of the magnetic field dependent trans- 
mission at fixed Fermi energy. However, in that case a 
semiclassical analysis of the quantum mechanical results 
is more involved for two reasons: Firstly, the areas of the 
trajectories are no longer unchanged by a reflection at 
the barrier. Increasing the barrier height leads to a shift 
of the spectrum to small areas |9|,[l7j . Secondly, in order 
to obtain well-resolved peaks in the area spectrum one 
has to perform the Fourier transformation over a rather 
large range of magnetic fields. Then, the condition that 
the cyclotron radius is much larger than the system size 
is no longer fulfilled in our numerical calculations and the 
shapes of the trajectories become field-dependent. This 
leads to broadening and thus to a severe restriction of res- 
olution. For high magnetic fields one may even find split- 
ting of the peaks in the area spectrum. Hence in most 
cases it will be hard to unambiguously identify peaks in 
the area spectrum. 

Low-frequency structure has been obtained in the 
(non-averaged) experimental area spectra of ballistic mi- 
crostructures [§]J|. The observed peaks correspond to 
areas close to those of the shortest periodic orbits. How- 
ever, a clear identification between peaks and trajectories 
has not been possible to establish. As discussed above 
the analysis becomes quite difficult when pairs instead of 
single trajectories are involved. Thus, a certain amount 
of caution has to be exerted for a detailed interpretation 
in terms of classical trajectories. 

We have shown that the length spectrum admits a sim- 
ple semiclassical analysis, even in the presence of tunnel- 
ing. In view of microwave experiments || as well as re- 
cent work on cavities in two-dimensional electron gases 
p2| which has demonstrated the possibility of measuring 
length spectra, it seems feasible to experimentally ver- 
ify the nonmonotonic dependence on the barrier height 
discussed here. 
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